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We report on form factors for the B — > Kl + l~ semi-leptonic decay process. We use several lat- 
tice spacings from a = 0.12 fm down to 0.06 fm and a variety of dynamical quark masses with 
2+1 flavors of asqtad quarks provided by the MILC Collaboration. These ensembles allow good 
control of the chiral and continuum extrapolations. The /b-quark is treated as a clover quark with 
the Fermilab interpretation. We update our results for /j| and f±, or, equivalently, /+ and /o- In 
addition, we present new results for the tensor form factor fj. Model independent results are 
obtained based upon the z-expansion. 
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1. Motivation and theoretical background 

The rare B meson decay B — > Kll is mediated by a flavor changing neutral current (FCNC). The 
Standard Model (SM) contribution occurs through a penguin diagram for the b — > sll process. Since 
FCNCs are small in the SM, this decay presents an opportunity to detect beyond-the-Standard- 
Model (BSM) physics. This decay has been studied by several experiments including BaBar [Jl|], 
Belle [g], CDF [§], and LHCb [Q]. We expect that LHCb will continue to improve its precision, 
and new B factories such as SuperB and SuperKEKB will greatly add to our understanding of this 
decay. 

Recently, the CDF Collaboration published results for the B + — > K + pL + pL~ differential branch- 
ing ratio [gp. They compared their measurement with a decade-old prediction of the decay form 
factors based on light cone sum rule (LCSR) calculations M\. There is considerable uncertainty in 
the LCSR form factors at small momentum transfer q 1 . Unless the uncertainty in the form factors 
is reduced, our ability to search for BSM physics will be severely limited by lack of knowledge of 
the SM decay rate. 

Lattice QCD enables first-principles calculations of the form factors with controlled errors that 
are systematically improvable. In fact, the first study of this decay, using the quenched approxima- 
tion, was done over 15 years ago [§]. Since then, there have been three additional studies using the 
quenched approximation [^, ^, g]. Recently, there have been two groups studying this decay using 
2 + 1 flavor improved staggered quark configurations from the MILC collaboration [0, [K 



2. Ensembles 

Table [T] shows the ensembles that were used in the current calculation. The MILC ensem- 



bles Qll] , m , |13| ] were generated using three flavors of dynamical quarks. The quarks are asqtad- 
improved staggered quarks and a Symanzik improved gauge action is used. The common mass (in 
lattice units) of the two light sea quarks is denoted ami and that of the strange sea quark is am s . 
We use four coarse (a « 0.12 fm), five fine ((a « 0.09 fm) and two super-fine (a « 0.06 fm) en- 
sembles, in order to control the chiral and continuum extrapolations. The ratio of light-to-strange 
sea-quark masses varies from 0.4 to 0.05. On each configuration, we use four source times to 
increase statistics. 

3. Form factors 

Semileptonic decays of a heavy-light pseudoscalar meson to a pseudoscalar or vector meson 
are characterized by form factors that describe matrix elements of a hadronic current between initial 
and final states. We concentrate here on the decay B — > Kll for which we need two matrix elements 
{K{pK)\isy tl b\B{pB)) and (K(pK)\isG >lv b\B(pB)), where the B meson momentum is ps and that of 
the kaon is pk ■ For the vector current matrix element, we define two form factors / + and fo, while 
for the tensor operator, we only need a single form factor fj : 

{K\isfb\B) = f + tf) U+ft-^L^kA + /0 (S 2 )^S^V, (3-1) 
V 1 / 1 
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f» a (fin) 


size 


ami /am s 


^rneas 


0.12 


20 3 x 64 


0.02/0.05 


2052 


0.12 


20 3 x 64 


0.01/0.05 


2259 


0.12 


20 3 x 64 


0.007/0.05 


2110 


0.12 


20 3 x 64 


0.005/0.05 


2099 


0.09 


28 3 x 96 


0.0124/0.031 


1996 


0.09 


28 3 x 96 


0.0062/0.031 


1931 


0.09 


32 3 x 96 


0.00465/0.031 


984 


0.09 


40 3 x 96 


0.0031/0.031 


1015 


0.09 


64 3 x 96 


0.00155/0.031 


791 


0.06 


48 3 x 144 


0.0036/0.018 


673 


0.06 


64 3 x 144 


0.0018/0.018 


827 



Table 1: Ensembles of asqtad Nf = 2 + 1 configurations analyzed. 



(K\iso» v b\B) = ^I^-y BP l- P l P ^ (3.2) 

In these equations, q^ = p^ — p K is the momentum transfer, and q 2 is the outgoing dilepton invari- 
ant mass squared. We can alternatively describe the vector current form factors by f\\ and f±. We 
define 

(K\isfb\B) = ^W B Yf\\ (Ek) + Pj_f±( E K)] , (3-3) 

where = /Mb is the four- velocity of the B meson and p±= p K — (pk ■ v)v^. 

In lattice QCD, it is convenient to work in the rest frame of the B meson. The new form factors 
are considered to be functions of the kaon energy Eg: 

, . (K\isfb\B) 
(K\isfb\B) 

a{Ek) ~ "VST (3 - 5) 

ME K ) = Ms^ K (K(mba-s\B(p)) 

where there is no sum on the index i in the equations for f± and fj. In the B-meson rest frame, 
<? 2 = {pb — Pk) 2 = Mg + M K — 2MbEk, so small Ek corresponds to large q 2 . 



4. Chiral-continuum extrapolation 

Because we work with light quark masses larger than in Nature and at non-zero lattice spacing, 
we must take the chiral and continuum limits of our results. We use heavy-light meson staggered 
chiral perturbation theory (HMS^PT) fll4) ] to perform extrapolations. We work to zeroth order in 
I /nib and next-to-leading order in the light-quark masses, kaon recoil energy, and lattice spacing. 
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Figure 1: The form factors /j| (left) and /j_ (right) as a function of Ek in lattice units. The legend in Fig. 2 
explains the meaning of all the symbols. 



The form factors f\\ and f± are fit to the following forms: 



f\\ 



C 



(0) 



fx 



1 + logs + C^mi + Cj z; (2m/ + m s ) + C^>E K + C^'E 2 + Cj 3 V 



C 



(0) 



f± = -h 



x_ 

f* 



1 



+ ■ 



logs 



-( 4 ) C-2 , n {5)„2 



Er+Ai+logs 



+ 



C 



(0) 



C^m, + (2m, + m s ) + C^ ; £jf + C^>E% + C£ 



^(2) 



(3), 



^( 4 )f2 , ~(5) 2 



(4.1) 



(4.2) 



where the notation in Ref. Q15| ] is used, and A# = — Mb- The 1-loop S?7 (2) chiral logarithms 
(denoted as "logs" above) are non-analytic functions of the pseudoscalar meson masses. We use 
the same form for fj as for f±, because the Isgur-Wise relation [ |l6| ] requires that fj = f + in the 
large Mb limit, and the shape of / + is dominated by f± at low recoil. Figures [I] and ||| show our 
results for fu , f± and fj including the continuum, physical quark-mass curve plotted as a black line 
with a cyan error band. 



5. z-expansion 



Because only relatively small kaon momenta are accessible in the lattice calculation, the range 
of q 2 directly calculable is close to the maximum value. However, the experimental results are best 
at low q 2 . We use the z-expansion to extend the kinematic range [17]. This expansion is based on 



the field theoretic principles analyticity and crossing symmetry. It is systematically improvable by 
adding more orders to the expansion. The q 2 region is mapped to z by 



yt+ — to 



\/t + -q 2 + y/U 



■to 



(5.1) 
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0.02/0.05 
0.01/0.05 
0.007/0.05 
0.005/0.05 
0.0124/0.031 
0.0062/0.031 
0.0047/0.031 
0.0031/0.031 
0.00155/0.031 
0.0036/0.0188 
0.0018/0.0188 
cont 




Figure 3: Fit of form factor multiplied by Blaschke factor and outer function to a polynomial expansion in 
z. Only statistical errors are shown. 

where t± = (Mb±Mk) 2 . We are allowed to choose to as we wish to minimize the range of |z|. It 
is convenient to choose to = t+ \ \ — yfl — j=J . This results in the full range of q 2 for this decay 
being mapped to \z\ < 0.16. The form factors are then expressed as a function of z: 

1 °° 

where the Blaschke factor B(z) = z{q 2 ,m 2 R ) is used to account for the pole structure of the form 
factor, and the outer function <j>(z) is selected such that Y,k=o a l — ^ t^l- 



In Fig. ||, we show our polynomial fits to the form factors multiplied by the Blaschke factor and 
outer function, i.e., f(q 2 )B(z)<p(z). To generate these plots, we calculate synthetic data at selected 
values of z in the range accessible to our simulations based on the chiral-continuum extrapolation 
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Figure 4: Preliminary results for form factors as a function of q 2 , including comparison with Ref. [gp. 

of the lattice data in Figs. [j]-|2| We see that the negative z region corresponds to maximum q 2 , 
and that is the region in which the lattice calculation is done. We impose the kinematic constraint 
f + (q 2 = 0) = fo(q 2 = 0) in the z-fit. We use terms up to z 3 for / + and fo, and z 2 for fj. The 
form factors as a function of q 2 are shown in Fig. ||. Both statistical and systematic errors are 
shown, but these results are preliminary and do not include the final renormalization factors or all 
the ensembles that we will analyze. 

At large q 2 , where we can directly compute the form factors, the error is about 5%. As q 2 
decreases, the errors from statistics, chiral fit and z-expansion grow to more than 30% for /+ and 
fo and to about 100% for fj. One reason that the relative errors grow at small q 2 is that the form 
factors themselves are getting much smaller. 



6. Conclusions 

We have presented some preliminary results for the rare decay B — > Kll form factors. There 
are additional ensembles not yet included in the analysis. We still need to include the perturbative 
renormalization factors, which are expected to be close to unity and will not change our results sig- 
nificantly. Our results for /j| , f± and fj can be used to calculate the SM B — > Kll branching fraction 
or to make predictions for this process in BSM theories. Other semileptonic decays currently being 
analyzed include B — > Ttlv and B s — >• Klv. 
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